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Abstract 

This paper presents a robust signal classification scheme for achieving comprehensive spectrum 
sensing of multiple coexisting wireless systems. It is built upon a group of feature-based signal detec- 
tion algorithms enhanced by the proposed dimension cancelation (DIC) method for mitigating the noise 
uncertainty problem. The classification scheme is implemented on our testbed consisting real-world 
wireless devices. The simulation and experimental performances agree with each other well and shows 
the effectiveness and robustness of the proposed scheme. 

1 Introduction 

The current research on spectrum sensing is mainly focused on detecting signals transmitted by licensed 
primary users (PU), such as TV broadcast, wireless microphone (WM), which is motivated by the regula- 
tor's stringent protection requirement. Due to the continuing innovations in new wireless technologies and 
business models, it is expected that more and more heterogeneous systems will share the same spectrum 
resources in the future. For example, apart from the primary TV broadcast, some other technologies, such 
as IEEE 802.22, ECMA-392, 3GPP LTE and cognitive PMSE JT) are targeted at exploiting the TV band 
white space (TVWS) for providing new services. As a result, the challenge arised is not only to protect 
the legacy PU, but also to optimize the coordination and coexistence among different secondary devices 
and networks, especially heterogeneous ones. In this context, in addition to the detection of PU's signal, it 
is also desired to use spectrum sensing for acquiring knowledge on other coexisting networks or devices. 
For this purpose, we propose a robust signal classification scheme which can identify the major primary 
and secondary wireless systems coexisting in TV band, which is based on the combination of a group of 
feature based signal detection algorithms using properly designed decision rules. 

We first review the reported signal classification schemes for spectrum sensing in literatures. In 
the authors focus on the classification based on signal's power spectrum density (PSD) shape. The cyclic 
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Figure 1: Signal classification based on parallel detection algorithms 



prefix (CP) is also taken into account for classifying OFDM signals which have similar PSD. Most of the 
targeted signals are actually not coexisting in the same frequency band, such as UMTS, DVB-S, HiperLAN 
and IS-95. In [[3]-[5j, classification schemes for OFDM signals are proposed which are based on the analysis 
of their cyclostationarity. Only the general signal models are discussed in these studies without considering 
the classification of standard-specific signals. In [6|, a classifier for LTE and WiMax is proposed based 
on detailed analysis of the cyclostationary features from both CP and pilot sub-carriers. In [7], network 
identification based on embedded cyclostationary signatures in OFDM signal is proposed. However, this 
scheme is not applicable to existing standards which don't have these signatures. 

The designed and implemented classifier presented in this paper covers the major wireless systems 
which exist or will potentially operate in UHF band below 800 MHz, they are DVB-T, 3GPP LTE, IEEE 
802.22, ECMA-392 and WM signal. We propose the dimension cancelation (DIC) scheme which can 
completely mitigate the noise uncertainty (NU) problem [[8j in receiver. Other harmful issues in sensing 
receiver such as unflat noise floor and spurs are also taken into account and mitigated. The proposed 
classifier is well validated with both computer simulation and experimental measurement using a testbed. 



2 The Signal Classification Scheme 

The proposed classification scheme is based on combining a group of feature based signal detection al- 
gorithms which can test the received signal in parallel (Fig {T]). Unique features are utilized by different 
detection algorithms, which makes the unmatched signal appears to be like noise. In this way, classifica- 
tion is achieved by matching the detection algorithms to received signal with properly designed decision 
rules. 
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2.1 Targeted Signals & According Detection Algorithms 

The targeted signals of the classifier are presented in Table [T] with their according detection algorithms and 
utilized features listed. It is well known that the goal of signal detection in spectrum sensing is to identify 
the following hypotheses: 

■Ho : y[n] = w[n] 

"H\ : y[n] = x[n] + w[n], 

where x[n] is the signal with channel's effect and w[n] is noise. For detecting the existence of the interested 
signal x[n], normally, a detection metric A should be defined, which is the result of certain processing on 
the received signal y[n]. The detection can be achieved by comparing the detection metric A with a 
threshold y which is predefined according to a desired probability of false alarm (PFA): 

A = T{y[n]} % y. (2) 

For detecting DVB-T signal, we adopt a robust detection algorithm called TDSC-MRC which was 
proposed in W\ and further studied by us in ] [T0| using real-world TV signal and with practical issues 
considered. It utilizes the autocorrelations based on the pilot subcarriers' period of four OFDM symbols: 

(N„-k)L-\ 

4=y= 2 y\n\y*[n + kL], (3) 



n=0 



in which y[n] is the received signal, L is the length of one period and N p is the number of periods within the 
observation time. Then the detection metric is composed by further correlating the adjacent autocorrelation 



values 



A 



TDSC-MRC 



Y.RiR, 



L 



k lv k+l 
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fc=i 

By analyzing the probability distribution of ^tdsc-mrc at *Hq with central limit theorem (CLT), the thresh- 
old according to the PFA P FA is derived as 



4 

yTDSC-MRC ~ 0~ w 



.1 



\ - Yi {N p - k){N p - k - l) lnPpA > (5) 



in which cr\, is the power of noise w[n]. 

For LTE signal, we adopt the algorithm utilizing CP based autocorrelation convoluted with a sliding 
windowed (SW) [11], which is named CP-SW detection here. First, the autocorrelation with time lag of 
OFDM's DFT size is performed: 

r[n] = y[n]y*[n + N DFT ]. (6) 
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Table 1: Signal Features and Detection Algorithms Used in Classification 



Standards/Types 


Algorithms 


Utilized Features 


DVB-T 


TDSC-MRC 9 10 


periodical pilot tone structure in time domain 


LTE (5 MHz mode) 


CP-SW (pre-align.) 11 


CP and DFT length, slot length 


IEEE 802.22 (8 MHz mode) 


CP-SUM 


121 


DFT length (indifferent to CP length) 


ECMA-392 (8 MHz mode) 


CP-SUM 


12 


DFT length (indifferent to CP length) 


WM 


PAR and DS of PSD 


high PAR & low DS in PSD due to narrowband 



The autocorrelation values r[n] are further aligned according to the OFDM symbol length 



I N-N DFT +l I , 



R[n] 



V r[n + l(N DFT + N C p)] n = 0, . . . , N DFT + N C p - 1. 



(7) 



/=o 



in which N is the total number of samples in the observation time. The detection metric of CP-SW is then 
formulated as 



Acp-sw - max 



i+Ncp-l 



2 R[n] 



i = 0,1,..., Ndft +N C p-l, 



(8) 



where R[n] is the cyclic extension of R[n] with CP length. If the signal has different CP lengths, such 
as the normal CP mode of LTE in which the first symbol of a slot has longer CP length than the other 
symbols, we propose to use pre-alignment on r[n] according to the periodical frame length: 



r'[n] = ^r[n + kL fnn ], 



(9) 



in which Lf m is the length of one periodical frame, for LTE it is the length of one slot. The probability 
distribution of A C p-sw in "Ho is difficult to obtain, hence we cannot derive the closed-form expression of 
the threshold jcp-sw- However, it can still be estimated through empirical method with large amount of 
detection tests in hypothesis "Tin- 

For ECMA-392 which use CSMA/CA MAC, the signal behaves as random bursts and has no peri- 
odical structure within the observation time. However, the detection metric can still be composed by the 
summation of r[n] |12j: 

AcP-5 UM 



N~N DFT -\ 

J r[n] 

n=0 



VN - Ndft + 1 

By analyzing the probability distribution of A C psum using CLT, we derive the threshold 

J cp-sum = cr 2 w V- In Pfa- 



(10) 



(11) 



For WM signal, we make use of its narrowband (< 200 KHz) feature exhibiting sparse distribution and 
high peak-to-average ratio (PAR) in estimated PSD of the signal in one TV channel. The PSD is estimated 
using the Welch's method: 



Y[m] = 



v £-i 
A-=o 



m = 0, 1,...,M- 1, 



(12) 



where M is the DFT size, D is the shifting step of DFT segment and the smoothing window v[n] of 
rectangular window is choose for achieving good frequency resolution. The PAR of the estimated PSD is 
calculated as 

<D = max(7[m])/Y, (13) 

in which Y denotes the mean value of the estimated PSD Y[m]. Since the targeted OFDM signals may 
also exhibit high PAR in their PSD if they experience frequency-selective fading, another property called 
the degree of sparsity (DS) is also used together with PAR to classify the WM signal: 

^ m\u. (14) 

y[m]>(l+p)Y 

The detection metric for WM signal becomes the logic result that if O and *F are both larger than their 
respective thresholds: 

A WM = (<D > cf> AND *F < t/0, (15) 

which is either or 1 . 



2.2 Decision Rules 

After the received signal are processed by the detection algorithms resulting detection metrics A, for 
OFDM class i and A WM for WM, the rules in (fT6l) are used to decide which signal is occupying the 
channel (^HofdmJc or ^wm) or the channel is vacant ("Kq). The WM signal is classified as long as the 



PAR and DS satisfies (15). Since the different OFDM detection algorithms has different threshold, we 
take Kilji as unified metrics with common threshold of 1 for finding the matched algorithm. Because 
each algorithms is specific to certain standards and modes, one signal can be classified as long as its 
corresponding algorithm is matched. 

'HoFDM,k '-k = argmax{Kilji), A WM = and max(A,/7 ; ) > 1 

i i 

"Hwm -Awm - 1 (16) 
*Hq : max(A ! 7y ! ) < 1 and Awm - 

i 

3 Mitigation of Receiver Imperfections 
3.1 Spurs and Unflat Noise Floor 

The narrowband spurs in practical receiver can be mistakenly regarded as WM signal by the classifier. 
Since the spurs' frequencies are normally fixed, the simple solution is to directly exclude the spurs' fre- 



quency components in the estimated PSD Y[m] when calculating PAR and DS values in ( 13 ) and ( 14 ). In 



the TDSC-MRC, CP-SW and CP-SUM detections which are based on time-domain autocorrelation, the 
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Figure 2: The structure of the spectrum sensing testbed 



spurs generate strong components in the autocorrelation which can degrade the detection and classification 
performances. The solution is to use band-stop filters to eliminate the spurs. 

The unflat noise floor can lead to nonuniform sensing performances for WM signals at different center 
frequencies. The simple solution is to pre-estimate the noise PSD W[m] at < H and use it to equalize the 
estimated signal's PSD: Y eq [m] = Y[m]/W[m] during the sensing tests. Then the Y eq [m] is used for further 
calculation of PAR and DS values. It is found that the autocorrelation based detection algorithms are not 
sensitive to the unflatness of noise floor. 



3.2 DIC for Mitigating Noise Uncertainty Problem 



It can be noticed in ^ and ([TT) that the thresholds for TDSC-MRC and CP-SUM aglorithms have the 



multiplication term cr^, or cr^. This means that the knowledge of noise power o\ is required in these 
algorithms. Practically, since cr^ cannot be exactly known due to the NU problem, the thresholds therefore 
become inaccurate, which leads to unpredicted PFA and probability of mis-detection (PMD). 
We found that using the new detection metric 

A' = A/(*5)° (17) 

instead of the original metric A, the knowledge of noise power a 2 w is not required. The & 2 = ^ 2Zn=o bWI 2 
is the estimated power of the received signal. The a is the factor which makes the denominator has the 
same dimension as the nominator in ( 17 ). For example, assuming the dimension of received samples y[n] 
is in volt(V), the detection metric A TDSC -mrc m <@> is V 4 , which requires a = 2. The metrics Acpsum m 
( 10) and A C p-sw in ^ have dimension of V 2 , hence, they require a = 1. This method is called dimension 
cancellation (DIC) in this paper. Since the PAR metric for WM signal is dimensionless, it is inherently 
unaffected by NU and the DIC is not needed. 



6 



Table 2: The Signals and Modes in Classification Test 



Standards/TVpes 


Modes 


DVB-T 


8 modes: all (he CP and DFT lengths 


TD-LTE DL (5 MHz mode) 


2 modes: TD-LTE (UL/DL conf. 2). long CP and normal CP 


802.22 DL (8 MHz mode) 


1 mode: 8 MHz, 1/16 CP 


ECMA-392 (8 MHz mode) 


1 mode: 8 MHz, 1/16 CP. duty cycle: 0.5 


WM 


1 mode: FM, loud speaker 



4 Simulation & Experimental Results 

The configuration of the spectrum sensing testbed is presented in Figj2} As listed in Table [2j five classes 
of signals with various modes are targeted in the tests. These signals are generated in Matlab according to 
their standard specifications emphasizing the features utilized by the classifier. They are also transmitted 
by E4438C in the sensing experiments. 

Three methods are employed for evaluating performances: 

• Method i(Ml): Use the signal and white noise generated in Matlab for simulation tests. 

• Method 2(M2): The signal is transmitted by E4438C and captured by USRP2+WBX with high 
SNR, the noise with spurs is captured by USRP2+WBX. Then they are combined according to the 
desired SNR in simulations. 

• Method 5(M3): Similar to M2, however, the received signal with local noise and spurs is directly 
used for sensing tests which features the real spectrum sensing scenario. The SNR is controlled 
using the ability of E4438C which can scale its transmitting power accurately. 

We first validate the effectiveness of applying DIC to TDSC-MRC, CP-SW and CP-SUM detections for 
mitigating NU problem with simulation. The NU is modeled using the "robust statistics" method [[8j, in 
which the upper limit of noise power is used to calculate the PFA while the lower limit is used to calculate 
the PMD. The simulation performances are presented in Figj3]which shows that the PFA performances of 
the three algorithms are notably degraded by NU of 1 dB. However, when DIC is applied, the detection 
performances become invulnerable to NU at all and are very close to the performances in the ideal case 
without DIC and NU. 

In the signal classification experiment, the sampling rate of the USRP2+WBX is set to 12.5 MS/s 
which can well accommodate one TV channel of 8 MHz. The signals are transmitted at center frequency 
of 560 MHz in which we've acquired transmission license from the local regulator. In the classification al- 
gorithm, we use observation time of 20 ms resulting 12.5M5 1 s x 20m^ = 0.25M samples being processed. 

For classifying WM signal, DFT length M = 256 and shifting step D = M/2 is chosen for estimating 
the PSD. We set the PAR threshold O = 2, the DS threshold W = 0.2 and p = 0.2 in ([13]) and based 



on empirical tests. For classifying the OFDM signals, the thresholds are obtained via both analytical 
calculation in ([5]) and ( [TT| ) as well as empirical tests with large amount of noise samples according to PFA 
of 0.01. The fixed sampling rate of 12.5 MS/s can introduce some out-of-band noise and signal which 
interfere the classification. Hence, before being processed by the detection algorithms for different OFDM 
standards, the received signal is filtered and resampled according to their standard- specific sampling rate. 
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Figure 3: Using DIC to mitigate NU for TDSC-MRC and CP based detections on DVB-T signal (2K, 1/4 
CP mode), observation time: 10 ms, PFA: 0.01 




SNR (dB) 



Figure 4: Overall classification performances of all the five signal classes (Table[2]) using Ml, M2, obser- 
vation time: 20 ms, PFA: 0.01 



The DIC is applied to all the OFDM detection algorithms. Besides, in M2 and M3 using the USRP2+WBX, 
the spurs at 560 MHz and 562.5 MHz are removed by excluding their PSD components in the WM detec- 
tion algorithm and by using band-stop filter in the OFDM detection algorithms. 

In performance evaluations, the five signal classes and their different modes are randomly chosen and 
processed by the classifier which gives the result of r H FDM,k, *Hwm or < K - The overall classification per- 
formance for all the signal classes are presented in Figj4} The probability of correct classification (PCC) 
characterizes how probable the existence of the signal is detected with its type correctly classified. Differ- 
ently, the probability of correct classification when detected (PCCD) characterizes the conditional proba- 
bility of correct classification when the signal is detected. In Figj4} the PCCD is close or equal to 1 at all 
the tested SNRs which shows a favorable property of the classifier that as long as the targeted signals are 
detected, they can be always correctly classified. 

FigJ5]further presents the PCC performances of the five individual signal classes. It shows good agree- 
ment between the simulation results using M1/M2 and the experimental results using M3 for all the signal 
classes, which well validate the effectiveness and feasibility of the proposed classification scheme in real- 
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Figure 5: Classification performances of the five individual signal classes (Table[2]) using Ml, M2 and M3, 
observation time: 20 ms, PFA: 0.01 

world spectrum sensing scenario. It also shows that the methods for mitigating spurs and unflat noise 
floor used in M2 and M3 cause little degradation to classification performance by comparing with the 
simulation results in Ml with ideal signal and white noise. 



5 Conclusion 

This paper presents a novel signal classification scheme configured for identifying the major existing 
and emerging wireless systems in UHF band below 800 MHz. The DIC method is proposed which can 
completely eliminate the influence of noise uncertainty problem. Apart from simulation, the classification 
algorithm is implemented on a spectrum sensing testbed considering practical issues of spurs and unflat 
noise floor. The experimental results agree well with the simulation results, which successfully validates 
the the effectiveness and feasibility of the proposed signal classification scheme. 
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